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During the directional solidification of peritectic alloys, two stable solid phases (parent and peri- 
tectic) grow competitively into a metastable liquid phase of larger impurity content than either 
solid phase. When the parent or both solid phases are morphologically unstable, i.e., for a small 
temperature gradient/growth rate ratio (G/v p ), one solid phase usually outgrows and covers the 
other phase, leading to a cellular-dendritic array structure closely analogous to the one formed dur- 
ing monophase solidification of a dilute binary alloy. In contrast, when G/v p is large enough for 
both phases to be morphologically stable, the formation of the microstructure becomes controlled 
by a subtle interplay between the nucleation and growth of the two solid phases. The structures 
that have been observed in this regime (in small samples where convection effects are suppressed) 
include alternate layers (bands) of the parent and peritectic phases perpendicular to the growth 
direction, which are formed by alternate nucleation and lateral spreading of one phase onto the 
other as proposed in a recent model [R. Trivedi, Metall. Mater. Trans. A 26, 1 (1995)], as well 
as partially filled bands (islands), where the peritectic phase does not fully cover the parent phase 
which grows continuously. We develop a phase-field model of peritectic solidification that incorpo- 
rates nucleation processes in order to explore the formation of these structures. Simulations of this 
model shed light on the morphology transition from islands to bands, the dynamics of spreading of 
the peritectic phase on the parent phase following nucleation, which turns out to be characterized 
by a remarkably constant acceleration, and the types of growth morphology that one might expect 
to observe in large samples under purely diffusive growth conditions. 



I. INTRODUCTION 

The spontaneous emergence of complex microstruc- 
tural patterns during the solidification of alloys is a sub- 
ject of both fundamental and applied interest Dur- 
ing directional solidification, a sample is pulled in an ex- 
ternally imposed temperature gradient G with a fixed 
pulling speed v p . This setup has been used extensively in 
fundamental studies of solidification patterns because it 
allows one to study their formation under well-controlled 
growth conditions. Depending on the type of alloy and 
the ratio G/v p , various patterns are possible. During 
monophase solidification of a dilute binary alloy, solute 
redistribution leads to a well-known morphological in- 
stability (Mullins-Sekerka instability below a critical 
ratio G/vp, and cellular or dendritic patterns are typ- 
ically formed. For nondilute alloy concentrations close 
to a eutectic point, two stable solid phases of different 
compositions can grow from a metastable liquid. In this 
case, the two phases cooperate and form lamellae or rods 
parallel to the growth direction (coupled growth). For 
off-cutectic compositions, coexistence between dendrites 
and coupled growth structures is also observed. 

Much less is known about microstructural pattern for- 
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FIG. 1. Schematic phase diagram of a peritectic alloy. C, 
concentration of impurity B; T m , melting point of pure A; T p , 
peritectic temperature. C p , C p /3, and C pa are the composi- 
tions of the liquid, f3 solid, and a solid that are in equilibrium 
at T p . ATn and AT^ are the nucleation undercoolings for 
a and [3 phases, respectively. Dashed lines are metastable 
extensions of the liquidus and solidus lines. 
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mation in peritectic growth M, despite the fact that 
many industrially important metallic alloy systems as 
well as ceramics such as the high-T c superconductor 
YBCO are peritectics. A schematic phase diagram of a 
peritectic AB-alloy (where B will be called the impurity 
for convenience) is shown in Fig. |l|. It contains a peritec- 
tic point, analogous to the eutectic point, at which two 
different solid phases, the parent (primary) and peritec- 
tic (secondary) phases, coexist with a liquid of higher 
composition than either solid phase. Above the peritec- 
tic temperature T p , the parent phase is stable and the 
peritectic phase is metastable, whereas below T p , the op- 
posite is true. For comparison, in a eutectic, both solid 
phases are stable below the eutectic temperature, and 
metastable above, and the impurity concentration in the 
liquid falls in between the concentrations of the two solid 
phases. For a sufficiently low G/v p ratio, a dendritic 
array structure of the parent or the peritectic phase is 
typically observed, and which of these two phases is se- 
lected depends on the alloy composition and G/v p ||. 
In contrast, for a high G/v p ratio morphological instabil- 
ity is suppressed. In this case, banded structures made 
up of alternating layers of primary and peritectic phases 
perpendicular to the growth direction are formed. These 
structures have by now been observed in various peri- 
tectic systems, including Sn-Cd {|J|, Sn-Sb |J, Zn-Cu 
||, Ag-Zn 0, and Pb-Bi ||,|. It is worth noting that 
eutecticlike coupled growth structures, which are quite 
distinct from banded structures, have recently been ob- 
served in the Fe-Ni system jlo| . Whether stable coupled 
growth is theoretically possible during peritectic growth 
has remained an open question for quite some time jjj, 
and we will address this issue elsewhere. Here, we focus 
primarily on banded structure formation and phenom- 
ena associated with the dynamical spreading of one solid 
phase onto the other. 

Recently, Trivedi has introduced a one-dimensional 
(ID) model to explain the formation of peritectic 
banded structures for purely diffusion-controlled growth 
The conceptual banding cycle assumed in this model is 
as follows. Consider a melt with homogeneous compo- 
sition Coo < C p being solidified starting from a flat (il- 
liquid interface in equilibrium. The rejection of impu- 
rities B into the liquid during solidification leads to the 
buildup of a solutal boundary layer. As a result, the 
interface temperature decreases, following the liquidus 
curve in the phase diagram. If is large enough, the 
interface temperature eventually falls sufficiently below 
T p for the peritectic phase to nucleate heterogeneously 
at the solid-liquid interface before the growth of the a 
phase has reached its steady state. The newly nucleated 
[3 phase rejects fewer impurities than the a phase. Con- 
sequently, the magnitude of the solutal boundary layer 
decreases and the interface temperature increases, fol- 
lowing now the /3-liquid coexistence line in the phase di- 
agram. If Coo is low enough, such that the corresponding 
interfacial temperature is sufficiently higher than T p , the 
a solid may renucleate again before the steady state is 



reached, and the cycle repeats. Therefore, this model 
predicts that bands can form only when the composition 
falls inside a narrow window in the hypoperitectic region 
(C pa < Coo < Cpp) whose width depends on the nucle- 

ation undercoolings AT^ and ATjy. 

The first attempts to validate this prediction experi- 
mentally yielded contradictory results. Directional solidi- 
fication experiments with Pb-Bi and Sn-Cd alloys seemed 
to show that bands also form in the hyperperitectic re- 
gion (C p p < Coo < C p ), in apparent contradiction with 
this prediction. An attempt was made to resolve this 
"composition range paradox" by incorporating convec- 
tion effects jl2| , assuming the existence of a fully mixed 
liquid of uniform composition outside a purely diffusive 
ID boundary layer of finite thickness. This model, how- 
ever, yielded a banding cycle and band spacings that are 
inconsistent with experimental results, hinting that this 
boundary-layer approximation (typically valid for strong 
convection) is inadequate to describe these experiments. 
Around the same time, careful serial sectioning of solidi- 
fied Pb-Bi and Sn-Cd alloys revealed that the seemingly 
banded structures are actually oscillatory treelike struc- 
tures connected in three dimensions fl3[| , and not discrete 
bands, thereby resolving experimentally this composition 
range paradox. Following this finding, a more accurate 
model was developed that assumes a planar solidification 
front, but incorporates a fully two-dimensional convec- 
tion flow field (lj]. This model successfully reproduced 
the observed oscillatory structures. 

Following these studies, experiments were conducted in 
thin tubes to reduce convection [0|. For tube diameters 
smaller than 1 mm, truly discrete bands indeed became 
observable inside a narrow composition range predicted 
by the ID diffusive growth model. Surprisingly, however, 
it was also observed that when the tube diameter was 
further reduced, "islands" of the (3 phase formed inside 
the matrix of the a phase, instead of discrete bands. This 
observation suggests that there is a microstructural tran- 
sition from bands to islands if the system size is reduced. 
It was also observed that islands tend to form more easily 
for initial compositions closer to C pa . In addition, some 
spatially chaotic patterns were observed in some exper- 
iments. The formation of these structures is controlled 
by a subtle interplay between the nucleation process and 
the competition between the growth of the nuclei and the 
preexisting phase. In this respect, the one-dimensional 
model may not always be adequate to describe this com- 
petition because it assumes an infinite spreading speed 
for the newly nucleated phase. Moreover, the 2D convec- 
tion model assumes a flat interface and is hence not well 
suited to simulate heterogeneous nucleation and spread- 
ing. In order to model accurately the formation of these 
different structures, a truly 2D model of interface evolu- 
tion is necessary. The particular difficulty of this problem 
is that the microstructure formation is controlled by an 
interplay between nucleation and growth of the differ- 
ent phases. No steady-state growth mode exists, which 
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makes the whole problem explicitly time-dependent. 

In this paper, we use a phase- field approach 15-2^1 
to investigate the formation of this class of banded mi- 
crostructures in a purely diffusive regime and a 2D ge- 
ometry. The phase-field method eliminates the need of 
explicit front tracking and thus greatly simplifies the task 
of numerically solving the equations of peritectic solidifi- 
cation that involve three-phase junctions. A phase-field 
model for peritectic growth has recently been proposed 
IE3I. Here, we use an alternative model that is closer to 



the eutectic model of Wheeler et al. 22 1 . 

We first investigate the spreading of the peritectic 
phase on the primary phase after a single nucleation 
event. We characterize in detail the dynamics of the 
three-phase junction during spreading and find a mor- 
phological transition from discrete bands of a and j3 
phases to isolated islands of (3 phase when the system 
size is decreased, in qualitative agreement with experi- 
ments. Moreover, our simulations enable us to under- 
stand physically the basic mechanism that underlies this 
transition. We then investigate the effect of multiple nu- 
cleations on microstructure formation in large systems 
by supplementing the phase-field equations with a phe- 
nomenological stochastic nucleation law. 

The remainder of this article is organized as follows. In 
Sec. II, we write down the sharp interface and phase-field 
models. Section III is devoted to the study of the equi- 
librium properties of the phase-field model and Sec. IV 
describes the simulation method. Results are presented 
in Sec. V, followed by a summary and conclusions in Sec. 
VI. 



II. MODEL 

A. Sharp-interface model 

The sharp-interface equations are given by 
d t C = D L W 2 C, 

v n (C L -C 1/ ) = -D L d n C L , 



T = T p + m v {C L -C p )-T v K 



(1) 



(2) 



(3) 



where C denotes the concentration of impurity B, and the 
subscript v labels the solid a and [3 phases. Equation (^) 
is the diffusion equation for the solute in the liquid with 
the solute diffusivity Dl- We have assumed that diffu- 
sion in the solid is negligible (one-sided model). Equa- 
tion (Q) expresses the mass conservation at the moving 
interface, with v n and d n denoting the normal velocity of 
the interface and the derivative normal to the interface, 
respectively. Finally, Eq. (||) is the Gibbs-Thomson con- 
dition at the solid- liquid interface, with K, m v , [i v , and 
r„ being the interface curvature, liquidus slope, kinetic 



coefficient, and Gibbs-Thomson constant of phase re- 
spectively. The Gibbs-Thomson constants Y v are defined 
by 



(4) 



where 7^ is the surface energy of the ^-liquid interface 
and L„ is the latent heat of fusion for phase v, both taken 
at the peritectic temperature. Young's condition 



laL^uL + 7/3it,3L + Jafitap = 0, 



(5) 



must be satisfied at the trijunction points where three 
phases meet, where t M „ is the unit vector parallel to the 
[i-v interface and pointing away from the trijunction. 



B. Phase-field model 

To distinguish between the three possible phases (liq- 
uid, a solid, and (3 solid), we follow a similar approach 
to that of Wheeler et al. for eutectic solidification by 
introducing two nonconserved order parameters (phase 
fields) <f) and tp. The first distinguishes between solid 
((/> = 1) and liquid (cp — —1), the second between the a 
solid (-0 = 1) and the (3 solid (tp — —1). The solid-liquid 
interface is defined by the level curve <p — 0, and the in- 
terface between the solid a and (3 phases is defined by 
the level curve tp — when </> is positive. One important 
difference from Ref. J22| is that in our model tp takes the 
well-defined value tp = in the liquid. This modifica- 
tion is necessary because, in the model of Wheeler et al., 
the equation of motion for ip becomes a simple diffusion 
equation in the liquid. This introduces an undesirable 
new time scale in dynamical simulations that is removed 
in the present approach. 

As a third dynamical variable we need the composi- 
tion C, which is a conserved field. We define the scaled 
composition 



s(r,t) = [C(r,t)-C pP ]/AC a 



(6) 



where AC„ = (C p — C pt/ ), v = a,/3, is the concentration 
jump at the ^-liquid interface at T p . 

In terms of these quantities, the equations of motion 
that govern the dynamics of the system are given by 



dtp 



SF 
5F_ 



dt 



M(0)V^ 



(7) 
(8) 
(9) 



where F is the dimensionless free energy of the system 
(i.e., the Hclmholtz free energy, divided by the product 
of the system size and a typical value of the free en- 
ergy density that sets the physical energy scale), M(4>) is 
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the mobility of the impurities, and t$ and are (fast) 
relaxation times for the phase fields. These equations 
are of the standard variational form known from out-of- 
cquilibrium thermodynamics. Note that, since SF/5c is 
the local chemical potential Eq. (|^) is simply the con- 
tinuity equation for the impurity concentration with the 
mass current J given by 



(10) 



If there are no fluxes across the boundary of the volume 
where F is defined, dF/dt < and Eqs. (@)-(||) imply 
that the dynamics drives the system toward a minimum 
of free energy. 

The free energy functional of the system is assumed to 
be of the form 



/(<^,c)}d 3 r. (11) 



Since F, (ft, and ip are dimensionless, the coefficients 
and Wif, have the dimension of length: they determine 
the width of the diffuse interfaces. The form of the free 
energy density is chosen such that there are two minima 
at tp = ±1 corresponding to the a (+) and (3 (—) phases 
for <p = +1. There is a single minimum in the liquid 
corresponding to cp — — 1 and ip = 0, and /(</>, f/', c) has a 
single minimum as a function of c for fixed values of (p and 
ip corresponding to the three equilibrium phases. A con- 
venient way to match these requirements is to construct 
a free energy density of the form 



A 



1 



f^, V, c) = -{c + AM) + -A 2 [l + h{4>)]h{ij)} 



\{BM<f) 
1 



;B 2 [l + ^)]/lW}+#) 



[l + h (<j>)]g(i>) + -[l 



ip 2 . (12) 



Here, A is a positive constant, and A\, A 2 , B\, and B 2 are 
functions of temperature. The function g is a double-well 
potential with minima at ±1, and the function h must 
satisfy h(±l) = ±1 and h'(±l) = in order to keep the 
minima of / at constant values of (p and ip, independent 
of the value of c. We take 



5 (0) = l/4-0 2 /2 + ^ 4 /4, 
h(4>) =3(0-0 3 /3)/2. 



(13) 
(14) 



The functions g(ip) and h(ip) are similarly defined. It 
follows trivially from Eq. ( JI4 ) that the bulk phase free 
energy densities are given by 



f L = f(-l,Q,c) = ^(c- A t ) 2 + XB 1: 



A 



(15) 



fot = /(1, 1, c) = -(c + A a + A 2 f - X(B 1 + B 2 ), (16) 



fp = f(l, -l,c) = ^(c + A 1 -A 2 



X(B 1 -B 2 ). (17) 



For the mobility function M(<fi), we take 



M(<j>) 



Dl 
2A 



(1-0). 



(18) 



With this choice, the diffusion coefficient of the impu- 
rity is a constant equal to in the liquid and zero in 
both solids, which corresponds to the so-called one-sided 
model. A standard asymptotic analysis of the sharp- 
interface limit of the present phase-field model pjtl shows 
that Eqs. (@)-(||) reduce as expected to Eqs. (0)-(||). The 
relation between the parameters in the two sets of equa- 
tions is given in the next section. 



III. PHASE DIAGRAM AND EQUILIBRIUM 
PROPERTIES 

By applying the well-known common tangent construc- 
tion to the bulk free energy densities given by Eqs. fll5|)- 
(|i7|), we can construct the equilibrium phase diagram 
of the phase-field model. The equilibrium compositions 
can be expressed in terms of the temperature-dependent 
functions Ai, A 2 , B\ and B 2 (see the appendix). How- 
ever, since there are only four functions, we can at most 
fit four lines out of six in the phase diagram (i.e., three 
pairs corresponding to a-liquid, /3-liquid, and a-(3 coex- 
istence). We may choose to construct, say, the two liq- 
uidus lines and the two solidus lines and leave the two 



solid-solid coexistence curves determined by Eq. ( A3) and 
Eq. (A4). Since we are interested only in the behavior 
of the system at temperatures close to T p , we assume for 
simplicity that the liquidus and solidus lines are straight, 
and that the concentration jumps at the solid-liquid in- 
terface are constant (liquidus and solidus are parallel). 
We can then choose A\ and A 2 as constants and B\ and 
B 2 as linear functions of the temperature. The corre- 
sponding expressions for the functions A\, A 2 , Bi, and 
B 2 are given in the appendix expressed in terms of the 
dimensionless temperature field 



T = 



(T - T p ) 
ImJAC'a 



(19) 



which is a measure of the temperature relative to T p nor- 
malized by the freezing range of the a phase. 

In the present model, there exists a temperature- 
dependent concentration c u such that, in the solid, the 
solid a (f3) phase is thermodynamically stable only if 
c < c u (c > Cu). By comparing Eq. ( [l6| ) and Eq. (|l7|), it is 
easy to show that c u is exactly midway between the two 
solid-solid coexistence lines. In order to avoid a phase 
transformation in the solid far behind the solid-liquid in- 
terface, we require c u to be independent of temperature. 
One way to achieve this is to make the solid-solid coexis- 
tence lines vertical by choosing suitable parameters. This 
difference from a real peritectic phase diagram is unlikely 
to change the qualitative behavior of the system. A phase 
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FIG. 2. Phase diagram for our model system. Dashed lines 
are metastable extensions of the liquidus and solidus lines. 



diagram for the model system used in our simulations is 
shown in Fig. |^. 

The equilibrium interface profiles connecting different 
phases can be obtained by solving the time-independent 
one-dimensional version of the equations of motion with 
suitable boundary conditions. Since the chemical poten- 
tial must be constant at equilibrium, the relation 



/i 



SF 



dc 



(20) 



can be used to eliminate the concentration field from 
Eqs. (0) and (g). The appropriate value of for a cer- 
tain temperature is obtained from the common tangent 
construction. The two resulting coupled ordinary differ- 
ential equations were solved numerically using a Ncwton- 
Raphson method on a one-dimensional grid of spacing 
Ax. For simplicity, we assumed = = W. Un- 
less otherwise stated, all the results below are obtained 
for Ax/W = 0.8, which provides a good compromise 
between computational efficiency and accuracy. The re- 
sulting equilibrium profiles, centered at the origin, for 
the phase fields and the concentration for a-L equilib- 
rium and /3-L equilibrium at T p are shown in Fig. For 
solid-solid equilibrium, the interface profile of tf> can be 
obtained analytically because (f> = +1 is a constant: 



ipo(x) = — tanh( 



V2W 



(21) 



In all cases, the concentration profiles are given by substi- 
tuting the equilibrium profiles (f>o(x) and ipo(x) obtained 
previously into Eq. (|20|). With the equilibrium interface 
profiles at hand, we can calculate the surface energies 
JaL, 7/3 Zo an d 7q/3, defined as the excess Gibbs free en- 
ergy per unit surface area. They are given by the expres- 
sions 



[W$(d x <p f + W$(d x ip ) 2 ]dx, (22) 
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FIG. 3. Equilibrium profiles for (a) a-L equilibrium and (b) 
/3-L equilibrium at T v (A = 2.5). The composition is scaled 
according to Eq. (ph. 



where (j>o and i/jq are the equilibrium profiles of the phase 
fields connecting phases fj, and v. The same formula for 
the surface energies can also be obtained by a matched 
asymptotic expansion p3] . For the solid-liquid interfaces, 
the surface energies are obtained by numerical integra- 
tion. For this purpose, it is more convenient to convert 
Eq. (||^) to a form without the gradients of the fields. 
Making use of the steady-state one-dimensional version 
of the equations of motion and the fact that /i is constant 
in equilibrium, we obtain after some algebra 



d_ 

tlx 



dx ^ 



A*c). (23) 



Now one can integrate Eq. (|23|) from — oo to an arbitrary 
x and make use of the expression for the equilibrium 
concentration profile and the bulk phase values to show 
that for the solid-liquid interfaces, 

\[Wl{d x cj>of + Wl{d x ^ Q f] 



-B 2 [i+KWHKWti]. 



(24) 



Here the upper and lower signs are for a-liquid and f3- 
liquid equilibrium, respectively, and A and B are defined 
in the appendix. For the solid-solid interface, the sur- 
face energy 7 Q ^ can be calculated exactly and is equal to 
2v / 2W v ,/3. 
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FIG. 5. Temperature range within which both a-liquid and 
/3-liquid equilibrium solutions exist, versus A. (Circles: limit 
of existence for a-liquid interface solution. Squares: limit of 
existence for /3-liquid solution.) 



Related to the surface energies are the two capillary 



lengths cLq and d$ defined as 



~ (A<v) a (0/i/0c) 



(25) 



which can also be expressed in terms of the Gibbs- 
Thomson constants T v by 



\m v \AC v 



(26) 



These are two of the physical length scales that are rel- 
evant in pattern formation in solidification problems. In 
real systems, the capillary lengths are microscopic and 
much smaller than all other physical length scales in the 
problem. Ideally, one would like to adjust the capil- 
lary lengths in the model to match the physical length 
scale ratios by choosing suitable model parameters. Since 
dfi/dc = A, it follows that d$ depends on A as 



A ' 



(27) 



Hence, to have small capillary lengths, one would like 
to increase A. However, A cannot be chosen arbitrarily 
large for two reasons. First, the surface tensions them- 
selves depend weakly on A for T ^ T p . As shown in 
Fig. ||, these variations amount to a few percent over the 
temperature range of interest when A is varied by a factor 
of 5. Secondly, the temperature range in which equilib- 
rium interface solutions exist also depends on A. More 
precisely, with a fixed value for A, the a-liquid equilib- 
rium solution does not exist if T is below a certain value, 
and the /3-liquid interface solution ceases to exist if T is 
above another value, because if T is too low or too high, 
the free energy density loses a minimum at ip equal to 
+1 or —1, respectively. We have estimated the range of 
temperatures in which both solutions exist for different 



A by finding equilibrium solutions at different temper- 
atures. The results are shown in Fig. |[ We can see 
that this temperature range becomes narrower when A 
increases. From now on, we fix A = 2.5 unless otherwise 
stated. This is a compromise between having a large A 
and and a sufficient working temperature range in which 
our two-dimensional simulations can be carried out. 

To check Young's condition, we performed two- 
dimensional simulations at T p on a square grid with 
Ax/W = 0.8. The equilibrium angles around a tri- 
junction were measured and found to be consistent with 
Eq. (||) to within a few degrees. 

For a moving interface, there are also nonequilibrium 
kinetic effects related to the attachment of atoms at the 
interface and solute trapping. Since we are mostly in- 
terested here in qualitative aspects of the growth mor- 
phologies, we have not analyzed all these effects in de- 
tail. We checked, however, by performing dynamical one- 
dimensional simulations that nonequilibrium effects only 
lead to a deviation from local equilibrium that does not 
exceed the Gibbs-Thomson effect caused by interface cur- 
vature in two-dimensional simulations. 



IV. SIMULATIONS 

For our simulations, we cast the equations of motion 
into a dimensionless form. For simplicity, we take = 
W^j, = W and = tv, = r. By defining the dimensionless 
variables 



r - t 



(28) 



and the new variable 



fL = c + A 1 h(<f>) + -A 2 [l + h(<t>)]h&), (29) 
the equations of motion can be written in the form 



G 



at 

dc 
dt 



of 



where 



^ ^ dip' 
aV • [£>(</>) V/i] 



tD l 

W 2 



(30) 
(31) 
(32) 

(33) 



TABLE I. List of simulation parameters. 



is the scaled diffusion coefficient of the impurity in the 
liquid, and 



= (1 - <j>)/2. 



(34) 



Instead of the concentration far from the interface, we 
may also use rip, the volume fraction of (3 formed in the 
solid, to characterize the overall composition of the sam- 
ple. The two quantities are related by 



Coo = (1 - rj )C pa + -qpCpp. 



(35) 



In a typical directional solidification experiment, the 
sample is pulled under a temperature gradient G with a 
pulling velocity v p . We define the dimensionless temper- 
ature gradient and velocity, G and v, by 



G 



G 



\m a \AC a 



(36) 
(37) 



Usually, thermal diffusion is orders of magnitude faster 
than the diffusion of the impurities, and hence we use the 
"frozen temperature approximation" which assumes that 
the temperature of the system adjusts instantaneously to 
the externally imposed temperature gradient. Accord- 
ingly, directional growth along the x axis is implemented 
by letting 



T = T + G(x — v p t), 



(38) 



where To is some reference temperature. 

There are five different physical length scales that con- 
trol the microstructural pattern formation: the two cap- 
illary lengths dp and dp defined by Eq. (gg), the two 
thermal lengths 
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and the diffusion length 

/ - Dl 



-W. 



(39) 



(40) 



Equations. (p0|)-(p2|) are integrated numerically on a 
two-dimensional grid. We use a = 1, Ax = 0.8, and 
At = 0.1. Zero-flux boundary conditions are applied to 
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the two sides that are parallel to the growth direction. 
There are several features in the model that can be ex- 
ploited in order to speed up the computation. First, the 
phase fields <p and ip differ significantly from ±1 only in 
the interfacial region, and hence we can avoid integrating 
Eqs. (|30| ) and (|!u]) away from the interface. In addition, 
Eq. Q32D needs to be integrated only in the liquid. Sec- 
ondly, the concentration field decays exponentially in the 
growth direction and varies only slowly in space in the liq- 
uid region far ahead of the interface. Hence, we can use a 
coarser and coarser grid as we move away from the inter- 
facial region. Thirdly, in order to simulate a semi-infinite 
system in the growth direction, we take advantage of the 
fact that all the fields remain unchanged in the solid in 
the one-sided model. Whenever the solid-liquid interface 
has advanced one lattice spacing, we pull the system back 
by one unit and keep the composition at the end of the 
liquid side at Cqq. With all these implementations, we 
are able to carry out simulations with typical lengths in 
the growth direction equal to about ten times the diffu- 
sion length. For the results presented in this article, we 
chose a pulling speed such as to have a diffusion length 
of Id — 200W. Other parameters and length scales are 
listed in Table §. 



V. RESULTS 

A. Dynamics of spreading 

Let us first concentrate on the spreading of the [3 phase 
on the a phase, starting from a single nucleus. Similarly 
to the situation considered in Trivedi's model, the simu- 
lation is started with a homogeneous composition in the 
liquid and a planar a-liquid interface. The lateral system 
size L is several times the diffusion length. Nucleations 
are assumed to occur heterogeneously at the solid-liquid 
interface when the temperature of the metastable inter- 
face reaches a certain undercooling with respect to the 
stable solid-liquid equilibrium. The nucleation under- 
coolings ATjy (a on (3) and AT|^ ((3 on a), shown in 
Fig. [I], are assumed to be constant. Accordingly, in our 
simulations a circular nucleus of (3 phase is put at the 
solid-liquid interface on one side of the box when the 
liquid composition at the interface reaches the threshold 
for nucleation fixed by the nucleation undercooling AT^. 
The radius of the nucleus is taken to be 6TU, slightly 
larger than the critical radius for nucleation. Since we 
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are interested here in the deterministic spreading dynam- 
ics following a single nucleation event, further nucleation 
is proh ibited . Multiple nucleation events will be treated 
in Sec. 



VD 



To characterize the dynamics of spreading, we recorded 
the position and velocity of the trijunction point. The 
sideways velocity v y can be regarded as a measure of 
the spreading speed of the (3 phase. Figures, ^(a) and 
||(b) show plots of Vy/vp versus time for different nucle- 
ation undercoolings and different compositions Cqq, re- 
spectively. Time is measured in terms of the diffusion 
time 



to 



D, 



(41) 



Two very different regimes of spreading can be clearly 
distinguished. Immediately after the nucleation, the 
spreading velocity is almost independent of the composi- 
tion, but strongly depends on the nucleation undercool- 
ing. The growth of the nucleus is influenced only by its 
immediate surroundings. On the length scale of the nu- 
cleus, which is much smaller than the diffusion length, 
the impurity concentration can be considered constant 
and is determined only by the nucleation undercooling. 
A higher AT^ is equivalent to a higher supersaturation, 
and hence a higher growth speed. 

At later times, the modifications of the diffusion field 
induced by the growing (3 phase influence the spreading 
dynamics, and the spreading velocities for equal under- 
cooling, but different Coo start to differ [Fig. §(b)]. Af- 
ter a complicated transient, the details of which depend 
on the choice of parameters, v y becomes a linear func- 
tion of time, which means a constant lateral acceleration 
of the trijunction. This acceleration is independent of 
the nucleation undercooling or the history of the system, 
but depends on the composition. An explanation of this 
finding can be deduced from Fig. [?], which shows the in- 
terface temperatures on the sides of the box and at the 
trijunction as functions of time. After the initial tran- 
sient, the temperature at the trijunction just follows the 
temperature on the a side. This implies that the a-liquid 
interface is almost planar up to the trijunction. We can 
also see that during the whole time of the simulation, 
the a-liquid interface is still relaxing toward its steady 
state below T p . Hence, the undercooling that drives the 
(3 phase to spread is increasing. 

In the late stages of spreading, the lateral diffusion 
length D/vy becomes much smaller than the solute 
boundary layer, and is comparable to or even smaller 
than the tip of the spreading finger. Therefore, the 
spreading speed should be a function of local supersatu- 
ration only. To check this assumption, we show in Fig. 
|(c) the same velocity curves as before, but now plotted 
against the undercooling of the planar a-liquid interface 
with respect to the peritectic temperature, —T. Since in 
our phase diagram the liquidus curves are straight lines, 
this undercooling is simply proportional to the super- 
saturation. The curves all collapse onto a single master 



curve after the initial transient, i.e., starting from the 
time when the interface ahead of the trijunction has be- 
come flat. This master curve is not linear, and does not 
smoothly extrapolate to zero. We did not attempt to 
calculate it theoretically. We expect that its detailed 
form should depend on the characteristics of the trijunc- 
tion, and in particular on the angles between the different 
interfaces. More theoretical and numerical work would 
be needed to elucidate in detail the role of the various 
material parameters. Remarkably, similar observations 
were very recently reported in experiments on a trans- 
parent organic eutectic alloy p7| during spreading of the 
secondary phase on a planar interface of primary phase. 
The spreading speed of the secondary phase showed an 
approximately linear increase with time, and the data 
could also reasonably well be rescaled onto an analogous 
master curve. 

Since the a-liquid interface far ahead of the trijunction 
stays fairly planar before the arrival of the (3 phase, the 
time dependence of the temperature on the a side can 
be well described by the Warren-Langer approximation 
p5[ . The rate of change of the supersaturation is solely 
determined by the composition Coo, which explains why 
the final slope of the curves in Figs, ^(a) and||(b) depends 
on Coo but not on the nucleation undercooling. 

A completely different behavior is observed when the 
composition is sufficiently low. As shown in Fig. ||(b) 
(dashed line), the initial spreading speed is the same as 
for the other runs. However, at later times, the spreading 
slows down and the trijunction point turns around such 
that v y becomes negative. Instead of a band, an isolated 
island of (3 phase is formed. This phenomenon will be 
addressed in detail in Sec. 



VB below. 



In the final regime of spreading, when the lateral speed 
becomes much larger than the pulling speed, the lat- 
eral diffusion length Dl/v v becomes comparable to the 
radii of curvature close to the trijunction point. In free 
growth, such conditions are reached only at very large so- 
lidification speeds. Under these circumstances, it is clear 
that the phase-field model does no longer reflect quan- 
titatively the sharp-interface equations, since it contains 
corrective terms due to the finite width of the interface. 
For instance, we observed a violation of Young's condi- 
tion at the trijunction point. More precisely, the angles 
between the interfaces, obtained by taking the tangent 
vectors to the = and ip = level curves at the tri- 
junction, are still consistent with local equilibrium, but 
the solid-solid interface is highly curved on a length scale 
comparable to the width W of the diffuse interface. This 
is due to the fact that the diffusivity varies smoothly 
within the diffuse interface, and hence the part of the 
solid-solid interface near to the trijunction is still able to 
move. As a result, the angles between the interfaces, seen 
on a macroscopic scale, differ from the local equilibrium 
angles. For the purpose of the present study, where we 
are mainly interested in the qualitative features of the 
microstructures, we did not investigate this effect quan- 
titatively. Let us remark, however, that such effects may 
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FIG. 7. Scaled temperatures at the trijunction and 
at the solid-liquid interface on both sides of the box. 
(AT^=0.03129, 770=0.375.) 



not be simply an artifact of the phase-field model, but 
may have a physical significance for high growth speeds 
if the relaxation of the trijunction toward local equilib- 
rium occurs on a time scale comparable to the time of 
diffusion through the trijunction region. 



B. Morphology transition 
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FIG. 6. Plots of Vy/vp (a) versus time for constant vol- 
ume fraction and varying nucleation undercooling, (b) ver- 
sus time for constant nucleation undercooling and vary- 
ing volume fraction, and (c) versus scaled undercooling 
— T of the aL-interface with respect to T p . Circles: 
77,3=0.375, AT£=0.01878; Squares: 77^=0.375, At£=0.03129; 
Diamonds: 770=0. 375, AT^=0.04381; Upward triangles: 
77^=0.25, AT^=0.03129; Leftward triangles: 770=0.3125, 
Af£=0. 03129. The dashed line in (b) is for 770=0.1875, 
ATjy=0. 03129, where an island is formed. 



The results of the preceding section were obtained for 
systems with lateral extensions of several times the dif- 
fusion length. For some sets of parameters, a surprising 
event occurs when the system size is reduced while all 
other parameters are kept constant. After some time, 
the spreading slows down, and the trijunction point may 
even turn around, such that v y becomes negative. As 
a result, the trijunction travels back to the wall where 
it originated, and an isolated island of (3 phase, or par- 
tial band, is formed. It hence appears that complete 
spreading is easier to achieve in larger systems, a quite 
counterintuitive result. Figsures ||(a) and ^|(b) show time 
series of typical snapshot pictures for the formation of an 
island and a band, respectively. The scales are the same 
on both axes and in both figures. Isoconcentration lines 
in the liquid are also shown. It can be seen that a lat- 
eral concentration gradient builds up in the liquid. This 
concentration gradient plays an important role in the in- 
terpretation of the morphological transition from islands 
to bands and will be discussed below. 

To study more systematically the conditions for the 
formation of islands, we performed simulations with var- 
ious lateral system sizes L and compositions Coo , with the 
following results 1 

1. At a fixed ATj^, there exists a critical composition 
c* such that if Coo < c* the /3 phase always forms 
islands. This critical composition decreases as the 
nucleation undercooling increases. 
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Spreading of a single j3 nucleus over the a phase 
= 0.02504, r)p = 0.3125. (a) Island formation at 
L/l D = 0.512, (b) Band formation at L/l D = 0.64. Time 
increases from top to bottom. (Dark region, a phase; light 
region, phase; unshaded region, liquid. The isoconcentra- 
tion lines in the liquid are evenly spaced in c.) 
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FIG. 9. Morphology map at (a) AT^ = 0.02504 and (b) 
Af£ = 0.03129. 



2. At a fixed ATj^ and if Coo > c*, there exists a 
critical lateral system size L c (c oc ) such that if the 
lateral system size L > L c the phase spreads 
completely and forms bands, whereas for L < L c 
it forms islands. L c decreases when either Coo or 
ATjy increases. 

3. When Coo < c* , such that the (3 phase always forms 
islands, the final shape (and also the size) of the 
islands is independent of L when L is larger than a 
certain size. 

Figure |^ shows the final morphology of the system 
for different Coo (or equivalently different rjp) and for 
different nucleation undercoolings. The dashed lines in 
Figs. ^|(a) and (b) represent an estimate for the critical 
system size L^Coo) for the transition from bands to is- 
lands. It can be seen that both L c and c* are smaller for 
higher Af£. 

The existence of the critical size L c can be understood 
by noticing that both a and phases have to reject impu- 
rities in order to grow, but the concentration jump at the 
a-liquid interface is larger than that at the /3-liquid in- 
terface. Since is the stable phase below T p , there exists 
a driving force for the phase to spread. On the other 
hand, as rejects fewer impurities than a, the impurity 
concentration in front of the phase rapidly decreases 
after the nucleation. This creates a lateral concentration 
gradient and hence an impurity flow from the a to the 
side as can be clearly seen in Fig. || This lateral im- 
purity backflow will accelerate the growth of a and slow 
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down the growth of (3 and hence there is a competition 
between the two phases. 

To be more precise, we can consider the following scal- 
ing argument. Let us assume for simplicity a constant 
spreading speed v s for the /3 phase. Then the time re- 
quired by the [3 phase to spread across the system is L/v s . 
On the other hand, impurities diffuse laterally through 
the system on a time scale of L 2 /D^. If L/v s < L 2 /Dl, 
the (3 phase is able to spread over the a phase before a 
significant impurity backflow can occur. If the opposite 
is true, the impurities have enough time to diffuse and 
the growth of the (3 phase is slowed down. Hence, the 
critical system size is given by 



Another way to interpret the above criterion is to note 
that the "diffusion speed" , which is roughly the speed of 
the impurities diffusing laterally through the system, is 
given by D^/L. If the diffusion speed is smaller than v s , 
spreading occurs, and Eq. ( |42| ) follows immediately. 

Clearly, the above argument is only qualitative. We 
have assumed a constant lateral spreading speed in 
Eq. ([42]) , although Fig. || shows that the spreading speed 
varies with time, and hence we can give no explicit ex- 
pression for v s as a function of composition and nucle- 
ation undercooling. However, we can see from Fig. ^jthat, 
for any given time, the instantaneous spreading speed in- 
creases with increasing nucleation undercooling and in- 
creasing volume fraction of (3 phase. This observation, 
together with Eq. ([l2|), allows us to understand qual- 
itatively the shape of the curves L c (coo) in Fig. |^. In 
addition, this criterion allows understanding of the strik- 
ing finding that spreading is easier in large systems that 
in small ones. 



C. Banding and island formation 

So far, we have concentrated on how a single (3 nucleus 
spreads on the a phase. It is natural to ask what happens 
if renucleation is allowed. This is a complicated problem 
since nucleation is an inherently stochastic phenomenon, 
which cannot be consistently treated within our deter- 
ministic model. However, we can try to gain some insight 
by incorporating nucleation phenomenologically. We will 
proceed in two steps. First, we treat repeated nucleation 
in small samples by deterministic rules to make contact 
with the recent experiments in the Sn-Cd alloy system 
Then, in the next section, we investigate the in- 
fluence of multiple stochastic nuclcations on the pattern 
formation dynamics in large systems. 

For solidification in small systems, it can be assumed 
that nucleation occurs predominantly at the container 
walls close to the solid-liquid interface. The density of 
nuclei and the nucleation rate are very rapidly varying 
functions of the composition. Therefore, it seems reason- 
able to assume that a nucleus will form as soon as the 



concentration in the liquid exceeds the threshold corre- 
sponding to the nucleation undercooling. Accordingly, 
we incorporate repeated nucleation by the following rules 
(the nucleation of a is handled like the nucleation of 
before, by placing a small circular nucleus at the solid- 
liquid interface). 

1. A nucleus of the new phase is placed at one side 
of the container as soon as the undercooling of the 
interface exceeds the nucleation undercooling. If 
both sides of the box reach the threshold at the 
same time, one side is chosen at random. 

2. Once a germ has nucleated, further nucleation of 
the same phase is prohibited until the germ has 
cither completely spread across the system or com- 
pleted the formation of an island. 

We consider now the two nucleation undercoolings as free 
parameters, and study different cases. When both AT^ 
and AT^ are large enough such that the newly nucleated 
phase spreads completely before the original phase is able 
to renucleate again, banded structures are obtained. For 
a smaller AT^, the f3 phase does not spread completely, 
but forms an island. The a phase overtakes the (3 phase 
and continues to grow until the nucleation threshold for (3 
is reached again. The islands of (3 phase form alternately 
on each side. This is a result of the history of the system: 
as a result of the formation of the previous island, the 
concentration of impurities is lower on the side where 
the last island occurred, and hence nucleation of [3 is 
favored at the other side. Examples of these banded and 
island structures are shown in Fig. [lO|(a) to Fig. [u](d). 
The scales on both axes in these figures are the same, 
but Fig. |lO|(e) has a different scale from Figs. |lO|(a)-(d). 
This last picture was obtained by a simulation at much 
smaller AT^, and the lateral size of the system is smaller. 
In this case, an oscillatory structure is obtained which 
tends to approach a coupled growth steady state after a 
complicated transient. 

These results are in good qualitative agreement with 
microstructures obtained in small samples of Sn-Cd alloy 
p3[ . In the experiments, islands tend to form always on 
the same side of the sample. We believe that this is due to 
a slight lateral temperature gradient across the sample, 
which is always present in experiments. 

D. Nucleation controlled microstructures in spatially 
extended systems 

Until now, we have mainly focused on the microstruc- 
tures formed in small samples. It is interesting to ask 
what kinds of structure are to be expected in large sam- 
ples in the absence of convection. This situation could be 
achieved either in quasi-two-dimensional thin samples, or 
in a microgravity environment. To model the microstruc- 
ture formation, larger scale computations were carried 
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out. However, in a spatially extended system, multiple 
nucleations are unavoidable and must be incorporated 
in a way that is consistent with the predictions of clas- 
sical nucleation theory. We chose to extend our model 
by incorporating the effects of multiple nucleation in a 
phenomenological manner. 

For the nucleation of the (3 phase on a planar a front 
(the same arguments also apply to the nucleation of a 
on (3), classical nucleation theory predicts the nucleation 
rate 



I = I e 



-AF*/k B T 



(43) 



where Iq is a constant prefactor (with dimension equal to 
the number of nucleations per unit volume per unit time) 
and AF* is the activation energy for heterogeneous nu- 
cleation. Assuming that the critical nucleus is a spherical 
cap on a planar substrate (the spherical cap model), AF* 
is given, respectively, in two and three dimensions by 



AF* 




0- (1/2) sin 20' 
167r(2+cosfl)(l-cost 
12 



2D 
3D, 



(44) 



where AFb is the difference between the bulk free energy 
of the (3 phase and of the liquid phase. The contact angle 
9, as shown schematically in Fig. [ll], is determined by the 
balance of surface tensions parallel to the substrate, 

laL = lap + J0L COS 9. (45) 

Assuming that the system is locally in thermodynamic 



equilibrium, it can be shown that AFb is proportional 
to (T — T p ) |26| , such that for a quasi-two-dimensional 
system the nucleation rate for f3 on a can be written as 



I = 



/ 2D exp[-A/(T-T p ) 2 ] if T<T p 
if T > T v 



(46) 



where A is a constant, and l%u has now the dimension of 
number of nucleations per unit time and per unit length 
of the interface. A similar expression with w = when 
T <T P holds for a nucleating on (3. The 3D form of AF* 
is used in deriving Eq. ( f46| ) since, in practice, the size of a 
nucleus is still much smaller than the thickness of a thin 
sample. Equation (^) determines the local nucleation 
rate and hence the probability per unit time of a nucleus 
forming as a function of the local temperature at the 
solid-liquid interface. 

Unfortunately, both experimental and theoretical esti- 
mates of the free energy barrier and the kinetic prefac- 
tor are scarce in the context of heterogeneous nucleation, 
since the actual values may depend on complicated de- 
tails of the interfacial structure. Since, in the present 
study, we focus on morphological aspects of the large 
scale structure, we decided to treat the two quantities 
as free parameters. Moreover, we want to compare the 
stochastic simulations to the deterministic runs of the 
preceding sections. Consequently, we may eliminate one 
of those two parameters by the requirement of recovering 
the rules used previously. That is, in the deterministic 
simulations a nucleus was put at the solid-liquid interface 
when it reached the predetermined nucleation undercool- 
ing. In the stochastic runs, nucleation should therefore 
occur with probability 1 for the same interface tempera- 
ture. This condition will lead to a relation between the 
prefactor and the energy barrier in the nucleation rate. 

To proceed, let us first specify how we treat nucleation 
in the simulation algorithm. The interface is scanned at 
a regular time interval At^ , and nucleation is attempted 
at points regularly spaced by a distance Asm along the 
interface. The nucleation rate may be rewritten as 



I = 



w(T) 
AtjyAsN 



where 
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w exp[-A/(T- 
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FIG. 12. Schematic plot of the dimensionless nucleation 
rate w(T) versus temperature for two different choices of wo 
and A. 



is a dimensionless function of the interface temperature. 
At each test point, a nucleus is generated with probability 
1 if w > 1, and with probability w otherwise. That is, if 
w < 1, a random number £ uniformly distributed between 
and 1 is drawn, and a nucleus is generated if £ < w. 
As before, the nucleus is spherical and has a size of 6W. 
A possible drawback of the procedure outlined above is 
that the actual nucleation rate depends on the values 
chosen for At^ and Asn- However, it is reasonable to 
assume that the microstructures should not depend too 
sensitively on the choice of these parameters as long as 
their values are much smaller than the time and length 
scales of the pattern formation process. 

Now we can relate the prefactor and the barrier in the 
nucleation rate. In the preceding sections, nuclei were 
introduced deterministically when the nucleation under- 
cooling was reached, that is, at a nucleation temperature 



-jv 
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N 



1 — mp/j 



(49) 



This implies that in the stochastic algorithm we must 
choose 



wpP) = wo exp[-A/(T£ - T p f] = 1 



(50) 



Given this condition, there is only one remaining free pa- 
rameter; we choose the dimensionless kinetic prefactor 
wo- This parameter controls the temperature range over 
which nucleations occur. Indeed, w(T) is a function that 
rapidly increases in a narrow temperature range around 
Tjf. Increasing wq and A simultaneously while respect- 
ing the constraint Eq. (j50|), the rise of the nucleation 
rates becomes sharper, as shown schematically in Fig. [l^. 
Now consider an cv-liquid interface during its transient, 
when the temperature at the interface is decreasing from 
above T p to its steady-state temperature. If the tem- 
perature range over which the nucleation rate increases 
significantly is narrow, all nucleation events will occur al- 
most at the same time when the interfacial temperature 




(a) (b) (c) 

FIG. 13. Microstructures obtained from simulations with 
stochastic nucleation events for L/Id = 2.048, r\p = 0.25, 
AT$ = AT& = 0.03129, and (a) w = 5 X 10 3 , (b) 
wo = 5 x 10 4 , (c) wo = 5 x 10 5 . Growth is from bottom 
to top; dark regions, a phase; light regions, f3 phase. 



coincides with Tjy. As a result, the mean separation be- 
tween the nuclei will be small. On the other hand, if w 
is small, nuclei appear with a broader spread in interfa- 
cial temperature, and the mean separation between the 
nuclei will be larger. Since the mean separation between 
the nuclei plays a similar role as the system size in a 
small sample experiment, we might expect a morphology 
transition from bands to islands as wo is increased. 

Figure |lj shows the microstructures obtained in simu- 
lations for small u>o, ranging from 5 x 10 3 to 5 x 10 5 . The 
lateral system size is about twice the diffusion length, 
and we use periodic boundary conditions in the direction 
perpendicular to the temperature gradient. The lateral 
spreading of multiple nuclei leads to a jagged morphology. 
Each V-shaped site in the figures indicates a nucleation 
event (either /3 on a or a on 0) . A transition from irreg- 
ular banded structures to islands can be observed as wo 
increases. Note, however, that nucleation events occur in 
bursts, leading to a spatial periodicity along the grow th 
direction that can be clearly distinguished in Fig. |l3](c). 
This means that there is still a "banding cycle," now 
consisting of layers of a two-phase composite structure 
(particulate structure) and layers of pure a matrix. 

The values for u>o in Fig. [13] are somewhat small. The 
attempt rate wq should be related to the rate Io in 

30 



Eq. (|43[) , which is typically about 10 nuclei/ciir s for 
heterogeneous nucleation in metallic systems 1TTJ] . Hence 
we investigated the microstructures formed with larger 
values of wq, ranging from 5 x 10 11 to 5 x 10 43 . In this 
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FIG. 14. Microstructures obtained from simulations with 
stochastic nucleation events: (a) wo = 5 x 10 1 , (b) 
wo — 5 x 10 26 , (c) wo = 5 x 10 43 ; other parameters as in 
Fig. |l3|. The unshaded region is the liquid. 



range, we always obtain island structures that look qual- 
itatively similar (Fig. [u]), and not too different from 
Fig. f3(c). This implies that the microstructures ob- 
tained with multiple nucleation events are not very sen- 
sitive to wq when vjq is larger than some critical value. 
Hence we would expect predominantly island structures 
in spatially extended systems. 

The last statement, however, is valid only for the quite 
restrictive assumptions made in our model. Most impor- 
tantly, we have assumed that the probability of nucle- 
ation depends only on the composition in the liquid, and 
not on the local geometry of the interface. This neglects 
the presence of grain boundaries and impurities (bubbles, 
inclusions) which can considerably enhance nucleation. 
Such heterogeneities broaden the distribution of the nu- 
cleation rate as a function of temperature, and would 
hence favor bands. The evolution of the grain structure 
could in principle be modeled by including the local crys- 
talline orientation as an additional order parameter. An 
interesting perspective is that the interplay between nu- 
cleation at grain boundaries and spreading might select a 
certain grain size, since for small grains a spreading phase 
can engulf and hence "heal" grain boundaries, whereas 
for large grains nucleation at the solid-liquid interface (as 
modeled in our simulations) may occur and lead to the 
formation of new grains. Such a study, however, is largely 
beyond the scope of this article. 



VI. SUMMARY AND CONCLUSIONS 

We have developed a phase-field model to investigate a 
class of banded microstructures that form during the di- 
rectional solidification of peritectic alloys under purely 
diffusive growth conditions. We focused on a regime 
of large thermal gradients and low pulling speeds where 
both phases are morphologically stable and the interface 
dynamics is controlled by a subtle interplay between the 
growth and nucleation of two competing solid phases, 
rather than by the morphological instability of one phase. 
We restricted our attention to a generic peritectic phase 
diagram that simplified both the models and the compu- 
tations, but our approach is in principle flexible enough 
to be extended to phase diagrams of specific materials. 

The two-dimensional simulations of this model have 
shed light on three main aspects of banding: the tran- 
sition from islands to bands that has been observed in 
narrow samples where convection is suppressed |i~3| , the 
associated dynamical spreading of one phase onto the 
other, and the type of structures that one would expect 
to form in wide samples under purely diffusive growth 
conditions that are not presently accessible in an earth- 
based laboratory, at least for the alloys investigated to 
date. 

We have shown that the transition from islands to 
bands can be understood in terms of a competition be- 
tween the lateral spreading of f3 on a and the diffusive 
backflow of rejected impurities from a to /3. This com- 
petition leads to the surprising result that bands tend to 
form more easily in wider samples, in qualitative agree- 
ment with recent experiments in the Sn-Cd alloy carried 
out in small samples in order to suppress convection [fj"3| . 
The critical system size L c at which this transition oc- 
curs depends on the nucleation undercooling for the j3 
phase that influences the spreading rate, and on the al- 
loy composition, with L c becoming infinite when the vol- 
ume fraction of the (3 phase falls below a minimum value 
necessary for band formation. The influence of other pa- 
rameters, in particular the form of the phase diagram 
and the equilibrium angles at the trijunction, has not 
been investigated in detail here. 

When the peritectic phase fully covers the parent 
phase, its spreading dynamics is characterized by a re- 
markably uniform acceleration of the moving trijunction 
that depends on the composition, but not on the nu- 
cleation undercooling. This acceleration originates from 
an increase with time of the local supersaturation (driv- 
ing force for spreading) associated with the relaxation 
of the planar parent phase ahead of the trijunction to its 
steady state below T p , together with a direct relationship 
between the instantaneous speed of the trijunction and 
this driving force that depends on the material proper- 
ties, but not on the history or the overall composition of 
the sample. Moreover, the relative angles between phase 
boundaries at the trijunction during rapid spreading de- 
part significantly from those prescribed by Young's con- 
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dition, indicating a strong departure from local equilib- 
rium. Both predictions might be experimentally testable 
in transparent organic eutectic systems that exhibit sim- 
ilar spreading transients before coupled growth is estab- 
lished @. 

The formation of multiple nuclei in wide samples (L 3> 
L c ) adds a stochastic element to the interface dynamics 
that renders the range of possible patterns even richer. 
One can nonetheless distinguish two basic types of struc- 
ture that can be understood within the framework of 
the single island to band transition in narrow samples 
[L ~ L c ). The first is a discrete banded structure made 
up of separate jagged bands that span the whole width 
of the sample. The second is a particulate banded struc- 
ture made up of approximate rows of particles (islands) 
of the peritectic phase embedded in the matrix of the par- 
ent phase. The banded (particulate) structure is natu- 
rally selected if the mean distance between nuclei is larger 
(smaller) than the critical sample width L c for the island- 
band transition, and, moreover, simulations reveal that 
the particulate structure is preferred if nucleation is as- 
sumed to follow a classical nucleation law. Even though 
we modeled patterns in wide samples with such a law, 
we expect the transition from a discrete to a particulate 
banded structure to be generally governed by the mean 
distance between nuclei even if other nucleation mecha- 
nisms (such as wall-induced nucleation and nucleations 
at grain boundaries) play a dominant role. Both types of 
structure could conceivably coexist in the same sample if 
nucleation conditions change during growth. 

There are a number of possible extensions of the 
present study. One is to investigate the patterns that 
form for a somewhat larger range of pulling speeds where 
the parent phase is morphologically unstable, but the 
peritectic phase is still linearly stable. Another is to in- 
corporate the influence of convection in a fully consistent 
way to make contact with experiments over a wider range 
of sample sizes, which is now possible within a phase-field 
context G8.E!1. 
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APPENDIX A: RESULTS OF THE COMMON 
TANGENT CONSTRUCTION 

The common tangent construction allows one to deter- 
mine the equilibrium composition for two-phase equilib- 
rium for given bulk free energies /„ of the two phases. 
For two-phase equilibrium, the bulk phases must have 
equal chemical potentials /i = df u /dc and grand poten- 
tials O = /„ — fic. Solving the resulting equations for our 
model bulk free energies, we find for solid-liquid equilib- 
rium 



(gi ± \B 2 ) 
" L ~ {At ± \A 3 ) 



- At =f A 2 . 



s ~ {At±\A 2 ) 

The upper (lower) sign is for the a (/?) phase. 
For solid-solid equilibrium, we get 



(Al) 
(A2) 



-ss 



Eh 
A 2 
Eh 

A- 2 



At - A 2 , 



At + A 2 



For convenience, we define 



and 



A=(At±-A 2 ) 



B=(Bt± -B 2 ) 



(A3) 
(A4) 

(A5) 
(A6) 



where the upper (lower) sign is for the a (/?) phase. 

In order to relate the parameters in our model to a 
physical system, let us write 



Bt — Bit + Bt 2 T, 
B 2 = B 2 t + B 22 T. 



(A7) 
(A8) 



Here, c and T are the scaled composition and tempera- 
ture, respectively, defined in the text. Let AC\, and m v 
be the concentration jump at the solid-liquid interface 
and the liquidus slope of phase v, respectively, at T p . 
Let r be the ratio ACp / AC a , then the parameters At, 
A 2 , Btt, Bt 2 , B 2 t, and B 22 are related to these quantities 
in the phase diagram by 



A x = \{l+r), 
A 2 = \{l-r), 
Btt = \(l + r) 



B 21 = -{l-r) 



(1 + r) 
(1 + r) 



(A9) 
(A10) 
(All) 

(A12) 
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TABLE II. Parameters for our model peritectic system. 





38.2 wt % 


A 1 


3.30745 x 10" 1 


Cpa> 


22.1 wt % 


A 2 


3.38509 x 10" 1 




33.0 wt % 


Bu 


-2.56790 x 10" 3 


m a 


-6.71865 K/wt % 


S21 


-2.62818 x 10~ 3 


mp 


-2.17 K/wt % 


S12 


-0.5 






S22 


0.0 



B12 — 
B22 = 



1 

4 

— I 1 - r 

2 V mp 



mp 
m a 



(A13) 
(A14) 



In order to have vertical solid-solid coexistence lines, we 
have to choose the parameters such as to make B22 van- 
ish. The parameters for our model peritectic system are 
listed in Table 0. 
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